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We consider the method of entanglement witness operator to verify genuine multipartite entan- 
glement for single-particle W states involving N parties. In particular, linear optical schemes using 
photo detectors and beam splitters are proposed to implement two different types of witness oper- 
ator in experiment. The first scheme that requires only a single measurement setting is shown to 
detect genuine multipartite entanglement for the overall efficiency beyond 1 — 1/N. On the other 
hand, the second scheme with N + l measurement settings achieves success at a significantly lowered 
efficiency than 1 — 1/N. 



PACS numbers: 03.67.Mn, 03.65.Ud, 42.50.Dv 



I. INTRODUCTION 

Entanglement, an element of central importance in 
quantum information processing, is still far from our 
complete understanding despite a great number of ef- 
forts for the past decades. Earlier works on bipartite 
entanglement [l[ have been naturally extended to multi- 
partite systems, and the multipartite entanglement has 
been identified as a resource to implement quantum in- 
formation processing to larger extent, such as error cor- 
rection [1 , secret sharing 0] , and one-way quantum com- 
puting 4], to name a few. Nevertheless, multipartite en- 
tanglement has been far less understood, and its gener- 
ation/verification in experiment seems more demanding 
than the bipartite ones [f|. 

One of the important issues regarding many-body 
quantum systems is to verify genuine multipartite entan- 
glement among parties in one way or another. Genuine 
multipartite entanglement is distinguished from other 
types of entanglement by the participation of all par- 
ties in quantum correlations, and it is particularly dis- 
tinct from biseparable entanglement: Consider an TV- 
party quantum system represented by the state of form 
|*)®|$), where |*) belongs to the Hilbert space iJ® m (m 
parties) and |$> to H® N - m (N - m parties). Although 
this system can show quantum correlations to some de- 
gree among m parties and N — m parties, respectively, 
it is far from true iV-party entanglement. In general, 
a mixed state p is called biseparable if it is a mixture 
of pure biseparable states, with different bipartitions al- 
lowed among component systems. Otherwise, the state is 
genuinely multipartite entangled. One example of bisep- 
arable states is the tripartite system represented in the 
number-state basis as 



Pl23= -\B){B\ X2 ® 

+\\B){B\x. 



\\B){B\ 23 



(1) 



where \B) is the Bell state of two parties and |0) the 
vacuum state of the third 61. 



Given an TV-body system, it is an important, but non- 
trivial, task to determine whether the system possesses 
genuine multipartite entanglement. One possible ap- 
proach to this problem is to set up the entanglement 
witness operator W @, H, H, EH H3 m sucn a wa V that 
the ensemble average trj/oVVl becomes positive or zero 
for all biseparable states yjj. If it takes a negative value, 
it thus becomes a clear signature of genuine multipartite 
entanglement. In this paper, we consider the detection 
of genuine multipartite entanglement for sing le-p article 
W states based on the entanglement witness [T3|. Par- 
ticularly, we are interested in the implementation of the 
idea in an optical experiment involving beam splitters 
and photo detectors. Note that the same problem was 
previously investigated, but in a limited context [3]. 
Specifically, Nha and Kim showed that the pairwise en- 
tanglement between arbitrarily chosen two modes can be 
detected regardless of the photo detector efficiency [14 1 
using the entanglement conditions derived for continu- 
ous variables [H, d, Ezl EH- Of course, it may indi- 
cate the multipartite entanglement structure of the W 
states to some extent [l9| . however, it is not a rigor- 
ous proof of genuine multipartite entanglement [6j. A 
counter-example is the biseparable state in Eq. (1), for 
which any two modes, when the third mode is traced out, 
have nonzero pairwise entanglement. 

In this paper, we consider two types of witness oper- 
ator for W states and propose how to implement those 
in practice for single-photon W states [20] - The first 
scheme, which requires only a single measurement setting 
without the need of a full state tomography, is shown to 
succeed in verifying genuine multipartite entanglement 
for the overall efficiency over 1 — 1/iV, where N is the 
total number of modes. On the other hand, the second 
scheme with N + l measurement settings is shown to 
achieve success even at a lower efficiency than 1 — 1/N. 

The proposed scheme in this paper is derived from the 
observation that the quantum fidelity of a given state p 
with respect to an entangled state \E) = U\P) is the 
same as the fidelity of the transformed state p' = U'pU 
with the product state \P). Here, the operation U is 
the entangling action on the product state \P) to obtain 
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the entangled state \E). Therefore, the fidelity measure- 
ment is alternatively achieved after the inverse unitary 
operation ft is performed on the state p. Note that 
a similar idea based on the time-reversed operations to 
demonstrate many-body quantum coherences has been 
well known in the nuclear magnetic resonance (NMR) 
community [2lj | . In particular, Lee and Khitrin recently 
used the time-reversed sequences of entangling opera- 
tions, i.e. disentangling operation, to verify the 12-spin 
" Schrodinger-cat" state in [2~3 ]. 

This paper is organized as follows. In Sec. II, we briefly 
introduce the entanglement witness operator and specify 
it for the class of W states. In Sec. Ill, a general idea 
of measuring the quantum fidelity, which is essential for 
the implementation of witness operator, is presented and 
applied to the case of single-photon W states. In Sec. IV, 
the proposed scheme involving only a single measurement 
setting is analyzed with some experimental inefficiencies 
incorporated, and in particular, the optimization of the 
detection scheme for the case of asymmetrical W states 
is discussed. In Sec. V, an improved scheme with N + 1 
measurement settings is presented along with a modified 
witness operator, and the critical efficiency for successful 
entanglement detection is shown to be significantly low- 
ered than 1 — 1 /N. Finally, the main results of this paper 
are summarized in Sec. VI. 



II. ENTANGLEMENT WITNESS OPERATORS 
FOR W-CLASS STATES 

In this section, we briefly introduce the entanglement 
witness operator to detect genuine multipartite entangle- 
ment and specify it for the class of iV-partite W states. 
Suppose that one has produced a certain multipartite 
state p, presumably entangled and most likely mixed due 
to experimental imperfections, close to a target pure en- 
tangled state \*S?). An witness operator W can then be 
constructed in a form 

W = al- (2) 

where the constant a is taken as the maximum possible 
overlap of a pure biseparable state \<f>) with the genuine 
multipartite entangled state \^>). Namely, 

a = max| 0)e Bl<0l*)| 2 , (3) 

where B represents the set of biseparable states [l(| ■ It is 
now straightforward to see trjpW} > for every bisep- 
arable state p = pes so that only genuine multipartite 
entangled states can take a negative value over the wit- 
ness operator W. This is generally true for mixed states 
due to the linearity of the witness operator and the con- 
vexity of separable states. 

Although it seems very demanding to find out the value 
of a for a given state |^) through optimization in Eq. (J3]), 
Bourennane et al. identified the constant a with the 
maximum Schmidt coefficient of with respect to all 



bipartite settings [10( . More concretely, in a fixed bipar- 
tition of all parties, say (fc-th bipartition), one may 
choose an orthonormal product basis \ij), where \i) and 
\j) belong to the Hilbert spaces of two parties, respec- 
tively. Then, the target state |^) is represented in the 
same basis as |\&) = EyCj^ and the maximum singu- 
lar value Amlx of the matrix Cjv can be evaluated for the 
fc-th bipartition. The same steps must be taken to obtain 
Amax for every possible bipartition and a corresponds to 
the maximum among those Amax- 

Let us now consider the class of AT-partite W states, 
which is represented by a form 

\W) = Cl |l, (),••• ,0) + c 2 |0, 1,- ■ ■ ,0) 

+ •■■ + cjv|0,0,--- ,1). (4) 

Then, the constant a in Eq. §5§ is evaluated as 

a= l-min{| Cl | 2 }, (5) 

where min{|ci| 2 } is the minimum among all \ci\ 2 . Given 
the number of parties N, the constant a takes the small- 
est value 1 — 1/N for the class of symmetric W states 
(q = l/VN for all i), which may be thus less demanding 
to detect than asymmetric W states. 

III. FIDELITY MEASUREMENT 

In this section, we first present a general idea to imple- 
ment the witness operator in experiment and then apply 
it to the class of single-photon W states. Once the con- 
stant a is identified in Eq. ©, the remaining task is 
to measure the fidelity between the reference state \$f) 
and the state p in question, namely, (^Ipl^) [23[, as 
the quantum average of the witness operator becomes 
tr{pW} = a- (*|p|*>. 

For the case of continuous variables, Kim et al. have 
proposed an experimental scheme to measure fidelity by 
mixing the two fields at a beam splitter and measuring 
the Wigner function of the output [24]. In particular, 
they considered homodyne detection to measure the fi- 
delity of two Gaussian field states (25|. However, this 
approach, which needs the preparation of two states to 
compare, may not be suited to the case of witness op- 
erator: The reference state \^f) in the witness operator 
is only a target and one shall have produced instead a 
mixed state p in reality. In other words, one has no ref- 
erence state to compare with the subject p. In this paper, 
another method to experimentally evaluate the fidelity, 
which does not rely on the preparation of reference state 
and proceeds only with the state p, is proposed as follows. 

In a number of cases, a many-body entangled state 
|\&) can be generated from an initial product state 
l^i) |<^2) • • ■ |0jv) which subsequently undergoes some en- 
tangling unitary operations collectively represented by f . 
Namely, 

|*)=f|0>EEf|0 1 )|0 2 )...|^). (6) 
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Then, the fidelity of our concern can be expressed as 

= WW), (7) 



where 



p = U^pU. 



(8) 



Consider the case that the initial product state is a col- 
lection of number states, i.e., \4>) — |^>i)|<fe) • ■ • \<j>N) — 
\n\)\n,2) ■ ■ ■ |tijv)< Given the multimode state p, the fi- 
delity is then reduced to the photon counting distribution 
of the transformed state p' = W pU. 

More precisely, the state p is first subjected to the 
inverse unitary operation f/T to create a new state 
p' = U< pU . Next, one measures the joint probability 
Pn 1 n 2 ---n N that the photo detector at mode i counts rii 
photons for the state p', that is, 



P, 



n 1 n 2 ---n N 



(ni\(n 2 \ ■ ■ ■ (n N \p'\m)\n 2 } ■ ■ ■ \n N ). (9) 



This probability corresponds to the fidelity of our con- 
cern. 

Let us apply the above idea to the case of single- 
photon W states. An arbitrary iV-partite TF-state of 
the form in Eq. can be prepared by injecting a sin- 
gle photon into an array of beam splitters as shown in 
Fig. 1 [nj. Namely \W) = LTjl } 0, - ♦ - ,0), where the 
unitary operator U is a series of beam splitter actions, 



U 



D 



{N-IN} 



■B 



{12}- 



The beam splitter operator 



Bijj + x\ transforms two adjacent modes {ay, ay+i} into 
Wp a 'j+i\ as 



l 3 + l 



sin 6j 
cos 0j 



- cost 
sin^. 



j 

a j+1 



(10) 



where cos 9i (sin 9i) denote the transmissivity (reflectiv- 
ity) of the beam splitter [26| . The coefficients Ci in Eq. ^ 
are then given by 

C\ = sin^i, 

C2 = cos6>i sin #2, 



cj = n? =1 cos^i sinOj, 



n 



JV-l 



(11) 



at the output. If a phase shift is additionally carried 
out with the amount <f>j at the j-th output mode, the 
coefficients become Cj = Cje~ l(l>j (j = 1, • • • , N), 

The above generation scheme suggests that the fidelity 
(W\p\W) = (1,0,- •• ,0|f7tp[/|l,0, ••• ,0) can be mea- 
sured by injecting a given state p to the series of beam 
splitter in reverse order to produce the state pU and 
then measuring the counting probability Pio - o that only 
the first detector counts one photon and the others no 
photons. 
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FIG. 1: (a) Experimental scheme for generating an arbitrary 
single-photon IF-state in Eq. Q. A single-photon at one in- 
put is injected into an array of beam splitters with the trans- 
missivity of the j-th beam splitter given by cos 9j . Then, the 
coefficients a in Eq. Q are given as those in Eq. (lll[l . (b) Ex- 
perimental scheme for detecting genuine multipartite entan- 
glement based on the IF-class witness operator. A given iV- 
partite state p, with each mode represented by bi, is inversely 
injected as shown and each photo detector (PD) counts the 
photons at the output. See main text. 



IV. DETECTION OF MIXED W STATES IN AN 
OPTIMIZED SETTING 

In this section, we analyze the detection of genuine 
multipartite entanglement of the VF-state in Eq. (d]) 
based on the scheme outlined in Sec. III. In particular, 
we take into account the experimental imperfections such 
as non-perfect single photon source and inefficient photo 
detectors. 

Suppose that one creates a single photon source with 
a success probability ps represented by a mixed state 
Single = Ps|l)(l| + (1 — Psr)|0><0| H3, which is injected 
to the beam-splitter array in Fig. 1. Then, the output 
state generated becomes 

P = p s \W)(W\ + (1 - Ps )\0 ■■ - 0)(0 ■■ - 0\, (12) 

and one wishes to detect the genuine multipartite entan- 
glement of this mixed state. 

When this output entangled state is subjected to the 
inverse array of beam splitters in the same configuration 
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as the one used in its generation (Fig. 1), the counting 
probability Pio-.-o becomes VPs, where r\ is the efficiency 
of the photo detectors. Now, the detection of genuine 
multipartite entanglement turns out to be a success under 
the condition 



a < rips, 



(13) 



where the constant a is specified in Eq. ([5]). In particular, 
for the symmetric W states, the above condition reads as 



1 



1 

N 



< VPs, 



(14) 



in which the requirement of the overall efficiency, VPS, 
becomes the least demanding. 

In the above analysis, we have considered the same 
configuration of the beam splitters to generate and de- 
tect the multipartite W states. That is, the pure W state 
component of the actual state p in Eq. (jT3J) is the same 
as the W state used as the reference in the witness op- 
erator W. On the other hand, for the case of the states 
p asymmetrical under permutations, one can try to op- 
timize the detection scheme by adjusting the reference 
state as follows. Given the actual state p in Eq. (TT2)) . 
which is a mixture of the target state \W) and the vac- 
uum state, one may choose a different witness operator 
W' with other W state \ W) in Eq. © as reference. Then 
the condition for the overall efficiency to detect genuine 
multipartite entanglement becomes 



|<W|W)I J 



< VPS, 



(15) 



where a' is the constant for the state \W) in Eq. ([3]). 
Now, given the target state \ W), it is desired to minimize 
the value rmTrepyp by finding out the optimal reference 
state \W). 

As an example, let us consider an asymmetric tripartite 
W state, \W a ) = 0,0) + ±|0,1,0) + 4^|0,0,1>. The 
constant a for this state is 3/4=0.75 from Eq. flS)), so 
if one uses the state \W a ) itself as the reference in the 
witness operator, the requirement becomes 0.75 < r]ps- 
On the other hand, by minimizing the value of ^ w " w ,^-i 

in Eq. (I5J), one finds 12 - %y/2 w 0.686 < r)p s with 
the symmetric W state \W S ) = 4= 1 1,0,0) + 4^|0, 1,0) + 

4^|0, 0, 1) as the reference, \W) = \W S ). 

For a fixed number of parties, N, however, one 
can readily show that the value of r^y^Typ is tightly 
bounded as 



a' 



1 



Kwiw)] 2 - 1 A ' 



(16) 



The equality in Eq. (TT6")) holds for the symmetric W 
states, \W) = \W) = ^|1, 0, • • • , 0) + ^|0, 1, • • • , 0) + 
■ • • + — ^=|0, • • • , 0, 1). Therefore, in the present scheme, 
the overall efficiency rjps must be greater than 1 — 1/A to 



detect genuine A-partite entanglement and the choice of 
symmetric W states is regarded as best from a practical 
point of view. 

As a final remark, let us address the problem of lo- 
cal phase shifts that may occur in the multipartite state 
under test. Suppose that the prepared state is in a form 



\W) = -4= (e-*|l, 0, • • ■ , 0) + e-^ |0, 1, • • • , 0) 
V A 



+ ... + e -^«| 0; ... ,0,1) ). 



(17) 



where fa is the local phase shift of z-th mode. If these 
phases are completely unknown, we cannot figure out a 
single experiment setup appropriately, as our scheme re- 
lies on the inverse operation to disentangle a given state. 
Even though the given state is genuinely multipartite- 
entangled, the measured fidelity |(W|W')| could be zero 
in the worst case. To resolve this issue, if it arises, one 
must place a phase-shifter at each mode in the detection 
scheme of Fig.l (b) before each mode enters the array 
of beam-splitters. By covering the whole range of phase- 
shift at each mode, one can maximize the fidelity, which 
will lead to a success in entanglement detection. As the 
phase shifts are locally performed, this does not affect 
our judgment on the entanglement structure of the given 
state. 



V. IMPROVED EXPERIMENTAL SCHEME 
WITH N + 1 MEASUREMENT SETTINGS 

In this section, we discuss a possible improvement of 
the experimental condition by considering a modified wit- 
ness operator. In Ref. [28[, Haffner et al. introduced the 
witness of the form 



W = al 2 - Q, 



where 



Q = \W N )(W N \ - \ B Si){BSi 



(18) 



(19) 



Note that the identity operator I2 in Eq. (fT8| refers to the 
Hilbert subspace corresponding to a total of two quanta. 
In Eq. QUI, the state \BS t ) = \0i)\W N -i) (i = 1, • • • , A) 
is biseparable where the i-th mode is in the vacuum state 
|0) and the other A — 1 modes in the W state. The mod- 
ified witness W' differs from the one in Eq. ^ by sub- 
tracting from the VF-state \Wn) the biseparable states 
\BSi), which all give the maximum overlap, 1 — 1/N, 
with |Wiv)- The constant a in Eq. (fT8| can be obtained 
by maximizing the expectation value of Q over bisepara- 
ble states as 



max 



I0)6B 



(faQ\fa. 



(20) 



The optimal biseparable state \ fa) — \a)\b) for the maxi- 
mum was identified as a form of \a) = cos 6\ |0, ■ • • , 0)^ + 
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FIG. 2: Critical efficiency e c 



for the modified 



l-(JV-l)/3 

scheme of Sec. V involving N+l measurement settings (cross), 
in comparison with the critical value e c = I — 1/N for the sin- 
gle measurement scheme of Sec. IV (circle), as a function of 
the number TV of modes. The genuine multipartite entangle- 
ment is successfully detected for the efficiency nps > e c in 
each case. The parameter (3 in the witness W' of Eq. (|18p 
was chosen as /3(N - 1) = 1 - 10" 3 . 



Let us now consider the detection of symmetric W 
states using the witness W'. With the same experimen- 
tal imperfections characterized by rj and ps as in Sec. IV, 
we obtain tr{p w W} = a - r)p s [l - (N - l)/3]. There- 
fore, the detection of genuine multipartite entanglement 
is a success for the efficiency ijps > 1 _(jy_ 1 ) |3 ■ As the 
constant a is given by the maximization procedure out- 
lined above for a fixed (3, one can finally obtain the im- 
proved experimental condition by minimizing the fraction 
e c = prn^zjTg over the parameter /?. We numerically 

checked that this fraction decreases as (N — 1)(3 — > 1 for 
a given number N. The result is plotted as a function 
of N with the value of (N — l)/3 = 1 — ICT 3 in Fig. 2. 
As clearly seen, the previous critical value of 1 — 1/N 
is significantly reduced in the modified scheme. For in- 
stance, when N = 3, the improved condition becomes 
e c w 0.515 < rip s for 2/3=1- ICT 3 and a ps 5.15 x ICT 4 . 



VI. SUMMARY 



smO^Wk) and |6) = cos6» 2 |0,- • • , 0) N - k + sm6 2 \W N ^ k ), 
respectively, for a fixed bipartition of {k, N — k} modes 
[28| . Therefore, a can be numerically evaluated by opti- 
mizing the expectation value over the parameters #i , 9 2 
and k for a fixed value of (3. 

To implement the modified witness operator W' in ex- 
periment, given a certain state p, we need to measure 
tv{pl 2 } and (BSi\p\BSi) (i = 1, • • • , N) in addition to 
(Wn\p\Wn)- The ensemble average ti{pl 2 } is simply the 
probability that the total quanta is at most two, which 
can be measured in the same experimental setup as Fig.l 
(b) with no extra efforts. In fact, it does not make a sig- 
nificant difference in our case if we take I (entire identity 
operator), instead of I 2 , for the witness W' in Eq. (fT8|) . 
On the other hand, the fidelity (BSi\p\BSi) can be mea- 
sured in a similar setup to the one in Fig.l (b) by injecting 
i-th mode directly to the photo detector and the other 
modes to the inverse beam-splitter array of the N — I 
mode W-stdXe |Wjv-i)- Therefore, a total of N + 1 mea- 
surement settings are required to implement the witness 
operator W'. 



In this paper, the approach based on the witness op- 
erator was considered to detect genuine multipartite en- 
tanglement for the single-photon iV-partite W states. In 
particular, the experimental schemes using the photo de- 
tectors and the beam splitters were proposed to imple- 
ment two different types of witness operator. In the first 
scheme involving only a single measurement setting, the 
threshold overall efficiency was found to be 1—1/N for the 
symmetrical mixed W states, which therefore becomes 
increasingly hard with the number of parties N. On the 
other hand, in the second scheme involving a total of 
N + l measurement settings, an improved condition was 
obtained with a significantly lowered critical efficiency. 
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